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We study ensembles of globally coupled, nonidentical phase oscillators subject to correlated noise, 
and we identity several important factors that cause noise and coupling to synchronize or desy- 
chronize a system. By introducing noise in various ways, we find a novel estimate for the onset of 
synchrony of the system in terms of the coupling strength, noise strength, and width of the fre- 
quency distribution of the natural oscillations. We also demonstrate that noise alone is sufficient to 
synchronize nonidentical oscillators. However, this synchrony depends on the first Fourier mode of 
a phase-sensitivity function, through which we introduce common noise into the system. We show 
that higher Fourier modes can cause desychronization due to clustering effects, and that this can 
reinforce clustering caused by different forms of coupling. Finally, we discuss the effects of noise 
on an ensemble in which antiferromagnetic coupling causes oscillators to form two clusters in the 
absence of noise. 
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Introduction. Synchronization describes the adjust- 
ment of rhythms of self-sustained oscillators due to their 
interaction Such collective behavior has important 
ramifications in myriad natural and laboratory systems — 
ranging from conservation and pathogen control in ecol- 
ogy [2| to applications throughout physics, chemistry, 
and engineering 0, 0] ■ 

Numerous studies have considered the effects of cou- 
pling on synchrony using model systems such as Ku- 
ramoto oscillators (5j. In a variety of real- world sys- 
tems, including sets of neurons [6] and ecological pop- 
ulations 0], it is also possible for synchronization to be 
induced by noise. In many such applications, one needs 
to distinguish between extrinsic noise common to all the 
oscillators (which we study in this paper) and intrinsic 
noise, which affects each oscillator separately. Conse- 
quently, studying oscillator synchrony can also give in- 
formation about the sources of system noise Q. Nakao 
et al. [9( recently developed a theoretical framework for 
noise-induced synchronization using phase reduction and 
averaging methods on an ensemble of uncoupled identi- 
cal oscillators. They demonstrated that noise alone is 
sufficient to synchronize a population of identical limit- 
cycle oscillators subject to independent noises, and sim- 
ilar ideas have now been applied to a variety of applica- 



tions 



The aforementioned papers characterize a system's 
synchrony predominantly by considering the probability 
distribution function (PDF) of phase differences between 
pairs of oscillators. This can give a good qualitative rep- 
resentation of ensemble dynamics, but it is unclear how 
to subsequently obtain quantitative measurements of ag- 
gregate synchrony [l3| . It is therefore desirable to devise 
new order parameters whose properties can be studied 



analytically (at least for model systems). 

Investigations of the combined effects of common noise 
and coupling have typically taken the form of studying a 
PDF for a pai r of coupled oscillators in a specific appli- 
cation 13, 14 1 . Recently, however, Nagai and Kori 15 1 
recently considered the effect of a common noise source 
in a large ensemble of globally coupled, nonidentical os- 
cillators, and they derived some analytical results as the 
number of oscillators N —¥ 00 by considering a nonlinear 
partial differential equation (PDE) describing the density 
of the oscillators and applying the Ott-Antonsen (OA) 
ansatz 16, 171. 



In this Letter, we consider the complex interaction be- 
tween noise and coupling. We first suppose that each 
oscillator's natural frequency (w) is drawn from a uni- 
modal distribution function. For concreteness, we choose 
a generalized Cauchy distribution = \ 7 2 + (J_ ulQ )2 > 

whose width is characterized by the parameter 7. (The 
case 7 = 1 yields the Cauchy-Lorentz distribution, and 
Wo is the mean frequency.) We then investigate the effects 
on synchrony of varying the distribution width. Taking 
the limit 7 — > yields the case of identical oscillators; 
by setting the coupling strength to 0, our setup makes 
it possible to answer the hitherto unsolved question of 
whether common noise alone is sufficient to synchronize 
nonidentical oscillators. 

We then consider noise introduced through a general 
phase-sensitivity function [24| , which we express in terms 
of Fourier series. When only the first Fourier mode is 
present, we obtain good agreement between theory and 
simulations. However, our method breaks down when 
higher Fourier modes dominate, as clustering effects [j| 
Il0j | imply that common noise can cause a decrease in our 
measure of synchrony. Nevertheless, we show that such 
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noise can reinforce clustering caused by different forms of 
coupling. Finally, we consider noise-induced synchrony in 
antiferromagnetic coupling, in which pairs of oscillators 
are negatively coupled to each other when they belong to 
different families. 

Globally coupled oscillators with common noise. We 
start by considering globally coupled phase oscillators 
subject to a common external force: 



dt 



K 



N 



i=i 



& - 9-) + o-Z(9 t ) P (t) , (1) 



where 9i and are (respectively) the phase and nat- 
ural frequency of the ith. oscillator, K > is the cou- 
pling strength, and p(t) is a common external force with 
strength a. As in Ref. jlEj . we will later assume that 
p(t) is Gaussian white noise, but we treat it as a general 
time-dependent function for now. 

The phase-sensitivity function Z(9) represents how the 
phase of each oscillator is changed by noise. Such a 
function can also be used for deterministic perturbations 
(e.g., periodic forcing). In the absence of coupling, one 
can envision that equation (JT|) is a phase-reduced descrip- 
tion of an A^-dimensional dynamical system exhibiting 
limit-cycle oscillations and which is then perturbed by 
extrinsic noise: 



^=F(X)+,G(X) P (t). 



(2) 



One can reduce (HJ) to a phase-oscillator system of the 
form f t = lo + ctZ(X) • G(X)^(*)> where Z ( X ) is thc 



phase resetting curve (PRC) [191 ] . In this case, 2(9) = 
Z(X)-G(X). 

We study the distribution of phases f(w,9,t) in the 
N — » oo limit. We assume jj>] that / evolves according 
to a nonlinear Fokker-Planck equation (FPE) : 

^ + ^ [{to + Kr Bin(9) + aZp(t)}f] = , (3) 

where r is the (complex) Kuramoto order parameter 
r(t) = J Q n exp(i9)f(io,9,t)d9doj. The magnitude 
|r| characterizes the degree of synchrony in the system, 
and the phase arg(r) gives the mean phase of the oscil- 
lators. To obtain an equation for r, we follow the ap- 
proach of Nagai and Kori [15( and use the OA ansatz 
that the phase distribution is of the form / = [1 + 
£~ =1 [(aexp(i6>)) n + (a*expH0))"], where a = a(w,i) 
is a complex- valued function. This form of /f r makes 
it possible to perform contour integration and obtain 
r(i) = a*{co - 7i,*). (See [3] for a discussion about 
multimodal /f r .) 

We express Z in terms of its Fourier series: 
Z{9) = ?f + Em=iKcos(m0) + b m sm(m9)} = c + 
Em.ilcm exp(im6*) + c* m exp(-im0)], where c m = (a m + 



ib m )/2. We substitute the series expansions into (J3]) to 
obtain 



dr K , 2N 

— = -jr + iuj r + —r (1 - \r\ ) 



+ inp ic r + cl+Yl r * n <+i + £ r " c »-i • ( 4 ) 

\ n=l n=2 ) 



To study the magnitude of r, we let r = v]4exp[i(wo* + 
</>)], where A and <fi are real. 

In the Supplemental Material (SM), we derive an FPE 
for the evolution of the PDF q(A,(f>,i). To obtain an 
expression for the time evolution of A, we integrate both 
sides of the FPE with respect to d> and perform averaging 
based on the assumption d, that the PDF Q(A,t) 
evolves slowly compared to the time scale of oscillations. 
The resulting FPE has a steady state given by 



Qoo(A) 



C 



exp 



[2h(A) +a 2 P 2 (A)}dA\ 



a 2 Pi(A) 



(5) 



where h(A) = (K - 2j)A - KA 2 , the function Pi (A) := 
Tlo9 2 A dt = |(1 ~ 4) 2 En=i^"(l c »l 2 )> the Unction 
P*{A) := ±fo9A%tdt + tSfg^dt = 1(1 - 
A )Y,n=i An ~ 1 i n ~ A )(\ c n\ 2 ), and C 1 is a constant ob- 
tained from the normalization Qoo = 1- 

Generalized Cauchy distribution of frequencies. We 
apply the above results to extend the theory developed 
in Ref. [HI to generalized Cauchy distributions of os- 
cillator frequencies. We set Z = sin(6*), so b\ = 1 
and all other Fourier coefficients vanish. This yields 
Px(A) = A(l - A) 2 /2 and P 2 (A) = (1 - A) 2 /2 0. 
The order parameter signifying the transition between 
synchrony and asynchrony adopted in Ref. (l5| is thc 
maximum of the PDF Qoq. To find where Qoo attains its 
maximum, we set Q'^ = 0. This yields 



2h(A) P 2 {A) P[(A) 



a 2 Pi{A) 2 Pi{A) 2 Pi{Af 



= 0. 



(6) 



Using our expressions for h(A), P\(A), and P 2 (A) then 



gives A — min < 0, 



so we need K+a /2 > 2j 



for synchrony. 

The aforementioned techniques can be applied to many 
scenarios. The case in which a = has been stud- 
ied @, and Ref. provides a detailed discussion for 
7=1. Let's consider the case K = in which uncou- 
pled, nonidentical oscillators are driven by noise. Several 
studies have considered a noise-driven ensemble of iden- 
tical oscillators 0, [13, 12 1 , but there has been much less 
work on nonidentical oscillators. We begin with the case 
Z{9) = sin(#) to simplify our expression for Qoo in equa- 
tion flU to obtain Qoo (A) = exp (- t^i-A) ) • We 
expect to observe a peak at A > for a 2 — 4j > and a 
peak at A = for a 2 — 47 < 0. We confirm this prediction 
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by simulating an ensemble of N = 10000 phase oscilla- 
tors evolving according to equation ([T]). We constructed 
the generalized Cauchy distribution for the natural fre- 
quencies using ujj = uj + 7tan 2^[2j — (N + 1)] for the 
jth oscillator [20(. In Fig. [TJ we compare the computed 
PDF Qoo{t) with histograms of A obtained from direct 
numerical simulations (i.e., stochastic simulations). Ob- 
serve that we obtain a peak at A > for a 2 — 47 ss 0.08 
but a peak at A = for a 2 — 47 » —0.04. We obtain good 
qualitative agreement between Qoo (t) and A, though the 
noisy nature of the system entails some mismatch be- 
tween theory and direct simulations. 

The increase in synchrony is gradual as a 2 — 47 changes 
signs. Accordingly, in addition to using the position 
of the peak to measure synchrony, we also use E(A) = 
J AQoo(A)dA. We show our results in the right panel of 
Fig. [TJ Using both theory and simulations, we see that 
E(^4) increases with the strength of the common noise 
and decreases with the width of the distribution. As the 
figure illustrates, even systems with only 50 oscillators 
already exhibit very good agreement for the expectation 
E(A). 




FIG. 1: (Left) Plots of the PDF Qoo(A). We calculate curves 
from the analytical expression ([6]), and we plot circles and 
squares are from 50-bin histograms of data obtained from one 
realization of direct numerical simulations. The solid curve 
and circles are for the case a = 0.4 and 7 = 0.02, and the 
dashed curve and squares are for a = 0.4 and 7 = 0.05. 
(Right) Plots of the measure of synchrony E(A) versus a. We 
obtain the curves from analytical calculations ([6]), and the 
circles and squares represent data from a temporal average of 
one realization. The solid curve and circles are for 7 = 0.02, 
and the dashed curve and squares are for 7 = 0.05. 

General phase- sensitivity functions. We wish to study 
the effects of noise via a general phase-sensitivity function 
2(9) rather than just 2(9) = sin(0). A sinusoidal phase- 
sensitivity function is overly simplistic but one can 
approximate many functions 2(9) using only a few terms 
in its Fourier series. Nakao et al. [9] showed for uncou- 
pled, identical limit-cycle oscillators that higher harmon- 
ics of 2 can cause oscillator ensembles to form clusters 
around a limit cycle and that increasing the strength of 
common noise makes the oscillators more sharply clus- 
tered (i.e., their phases reside in a smaller interval). 
Equally-spaced (or almost equally-spaced) clusters lead 
to cancellation effects and a decrease in the value of the 
order parameter |r(i)|, which is problematic for our previ- 
ous analysis. Moreover, the formation of multiple clusters 



causes the OA ansatz to break down: from the normal- 
ization of the phase distribution /(w, 9, t), we know that 
\a\ < 1, so the coefficients of higher modes must have 
smaller magnitude than that of the first mode. This re- 
sult does not correspond to a phase distribution with 
multiple clusters. 

To illustrate the breakdown of the OA ansatz, we con- 
sider the example 2(9) = sin(2#) + cos(20). This func- 
tion can arise in an ensemble of Stuart-Landau oscilla- 
tors from adding multiplicative common noise (where the 
noise strength is multiplied by a function of one or more 
system variables). This yields Pi = A 2 (l — A) 2 /2 and 
P 2 = A(l — A) (2 — A)/2, which we insert into equation 
([5]) to calculate the steady-state pdf Q 00(A) and the or- 
der parameter E(^4). We also estimate the level of syn- 
chrony in the absence of noise by setting a = 0. (We use 
the notation Aq to denote values of A in this situation.) 
This yields H(Aq)Q = constant, implying that Q(Ao) di- 
verges at the zeros of h(Ao), which occur at Aq = and 
Aq = 1 — 2j/K. We show our numerical results in the left 
panel of Fig. O Observe that the presence of the higher 
harmonic leads to a decrease in synchrony rather than 
an increase in synchrony with increased noise strength, 
in contrast to many studies of noise-induced synchrony 

To characterize this decrease in synchrony, we use 
a family of order parameters from Ref. [HI to study 
clustering. We define A 2 — \r 2 1 2 , where r 2 (t) — 
J™ Jq* exp(2i6)f(u},8,t)d6du}. (One can similarly de- 
fine A q for all q € Z.) For the OA ansatz to hold, one 
needs A 2 = Aq. As we show in the right panel of Fig. [3] 
using direct numerical simulations, we find a high correla- 
tion between the clustering effect quantified by (A 2 
and the noise-induced decrease in synchrony quantified 
by Aq — (Ai). (The notation (x) refers to the temporal 
average of the variable x.) 
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FIG. 2: (Left) The synchrony measure ¥,(A) versus a for 
K = 0.5 and 7 = 0.1. The dashed curve is our erroneous 
calculation of synchrony, the solid curve is our estimate of syn- 
chrony in the absence of noise, and the circles are from direct 
numerical simulations. (Right) Comparison of the clustering 
effect (^2) — Aq and the noise-induced decrease in synchrony 
Ao — (Ai) from the left panel. 

Clustering. We now show that noise increases clus- 
ter synchrony when there is higher-order coupling (i.e., 
when the dominant mode in the coupling function is not 
the q = 1 Fourier mode). We take 2 q = a q sm(q9) + 
b q cos(q9) = c q exp {qi(9}+c* exp {— qi(9)} (where geZ) 



4 



to obtain 
dt 



K 
N 



TV 



^smiqlOj - 6$ + aZ q (9)p(t) , (7) 



which is discussed for the case a = 
Ref. M- 



in 



By defining the mode-g order param- 
eter r q (t) = f J Q 2Tr exp(qid)f(u},9,t)d8duj, we de- 
rive a similar nonlinear PDE as before: ^ + 
■§ s {[u} + Kr q sin(q0)+aZ q p(t)]f} = 0. Applying the 
same method as before yields 



dr q 
dt 



K , 



\r q \ 2 ) + iap{c* + r 2 q c q ) 
(8)' 

Setting r q = A q exp(qi6) and following the previously dis- 
cussed procedure yields the steady-state PDF Qoo(A q ) = 

{l mA t^ )]dA ' )^ where P l{ A q ) = 



7W 6XP 

A q {l - A q ) 2 {\c q \ 2 )/2 and P 2 (A q ) = (1 - A q ) 2 (\c q \ 2 )/2. 
This, in turn, implies that noise and coupling both in- 
crease the "g-cluster synchrony" of the system. We verify 
this for two and three clusters in Fig. [3J 
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FIG. 3: Cluster synchrony induced by noise and coupling 
via the phase-sensitivity function Z — sin(q9) + cos(q9) and 
coupling between oscillators of the form sin(q[8j — 8i]). We 
use the parameter values K = 0.1 and 7 = 0.05. The solid 
curves are from analytical calculations, and the circles are 
from direct numerical simulations. The insets show snapshots 
of 500 oscillators for a = 0.8. The left panel is for q = 2, and 
the right panel is for q — 3. 

Antiferromagnetic coupling. We now consider inter- 
actions of noise and coupling for antiferromagnetic cou- 
pling, in which there are two groups of oscillators with 
positive coupling between oscillators in the same group 
but negative coupling between oscillators in different 
groups. We label the two groups as "odd" and "even" 
oscillators. The temporal evolution of the phase of the 
ith oscillator is 



d(k 
dt 



1 N 
< 1 -^>" Si 



l )+aZ(0 l )p(t), (9) 



where Kij = K if i + j is even and = —K if it is 
odd. If a = 0, then the oscillators form two distinct 
clusters when K > K c = 27 (see the SM for the deriva- 
tion). If one defines an antiferromagnetic order param- 
eter r a f(t) = (1/N) J2j(~ 1) J exp(i#j), then the depen- 
dence of |r a f| on K and 7 is analogous to what occurs 



in the conventional Kuramoto model. However, as we 
have seen previously, the effect of noise when the first 
Fourier mode of Z dominates is to synchronize the os- 
cillators (i.e. to form a single cluster). We explore this 
using direct numerical simulations in Fig. 2J 

In agreement with our intuition, the noise and coupling 
have contrasting effects. Accordingly, |r a f| decreases with 
increasing o (see Fig. @h.). The noise causes the system 
to "jump" between states with low and high values of 
|r a f|, as shown in the inset. (By contrast, as shown in 
Fig. 0b, |r a f I decreases with increasing natural frequency 
distribution width parameter 7.) Furthermore, the de- 
crease in synchrony |r a f|o. = o ~ l r af| correlates positively 
with the increase in the traditional measure of synchrony 
M = 1(1/^0 Ej ex P(^) I = v^4 (seeFig.Sb). (The Pear- 
son correlation coefficient between |r a f | CT= o — \ r af\ and |r| 2 
is 0.955.) There is no such relationship in the case where 
7 is increased, and \r\ remains small and approximately 
constant (see Fig. [4}i). 
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FIG. 4: Results of direct numerical simulations for antifer- 
romagnetically coupled phase oscillators, (a) Antiferromag- 
netic synchrony |r a f| versus noise strength a for K = 0.05 
and 7 = 0.008. In the inset, we show a sample realization for 
a = 0.5 between times t = 1000 and t = 2000. (b) Antiferro- 
magnetic synchrony |r a f| versus 7 for K — 0.05 and a = 0.01. 
In the inset, we show a sample realization for 7 = 0.05 be- 
tween times t — 1000 and t — 2000. (c) Circles give the 
decrease of antiferromagnetic synchrony |raf| CT =o — kaf|, and 
crosses give the square of the usual Kuramoto measure of 
synchrony \r\ 2 . (d) Same as panel (c), except the horizontal 
axis is the natural frequency distribution width parameter 7 
rather than a. [Each data point in the figures in the main 
panels represents the temporal average of ((9]l with TV = 500 
oscillators.] 

Conclusion. We have examined noise-induced synchro- 
nization, desynchronization, and clustering in globally 
coupled, nonidentical oscillators. We demonstrated that 
noise alone is sufficient to synchronize nonidentical os- 
cillators. However, the extent to which common noise 
induces synchronization depends on the magnitude of 
the coefficient of the first Fourier mode. In particular, 
the domination of higher Fourier modes can disrupt syn- 
chrony by causing clustering. We then considered higher- 
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order coupling and showed that the cluster synchrony 
generated by such coupling is reinforced by noise if the 
phase-sensitivity function consists of Fourier modes of 
the same order. 

One obvious avenue for future work is to develop a the- 
oretical framework that would make it possible to con- 
sider multiple harmonics of both the coupling and the 
phase-sensitivity functions. It would also be interest- 
ing to consider generalizations of antiferromagnetic cou- 
pling, such as the variant studied in Ref. [2l| . One could 
also examine the case of uncorrelated noise, which has 
been studied extensively 22J via an FPE of the form 
^ + JL [oj + Krsm(9)} = However, proceeding us- 
ing Fourier expansions like the ones discussed in this Let- 
ter could perhaps yield a good estimate of the effect of un- 
correlated noise on such systems. Because of the second 
derivative in this system, the OA ansatz no longer ap- 
plies, and a generalized or alternative theoretical frame- 
work needs to be developed. 
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SUPPLEMENTAL MATERIAL 
Derivation of the Steady-State Probability Density Function for Q(A) 

Using equation we obtained an equation describing the temporal evolution of the complex order parameter r. 
To study |r|, we let r = \/AexTp[i(uJot + (/))], where A and <j> are real. We express the Fourier coeffiicients of Z in terms 
of their real and imaginary parts using c m = (a m — ib m )/2 and then collect real and imaginary terms to get 

dA 

— = h{A)+ag A {A,u o t + 0) Pl (10) 
d(f> 

— = ag ctj (A,u} t + <j))p, (11) 

where 

h(A) = {K- 2 7 )A - KA 2 , 

OO 

g A (A, uj t + (f>) = Yl An/2 { 1 - A ) {«« sin(n[cj i + </>]) - b n cos(n[cj t + <£])}, 

n=l 

j X 

g^A, u} t + 4>) = -a + A {n - 2)/2 {\ + A) {a n cos(n[uj t + 4>]) + b n sin(n[w i + <f>])} . 

n=l 

Thus far, we have not made any assumptions about the form of the external driving function p(t), but we now set it 
to be Gaussian white noise. If the correlation times of the noise is comparable to the amplitude relaxation time of a 
limit-cycle oscillator, then one might need additional terms to describe the exact phase dynamics [23| . However, such 



terms do not affect long-time phase diffusion and synchronization ll| . 

Treating equations (fTQ|) and (fTT|) as Ito stochastic differential equations (SDEs) yields a Fokker-Planck equation 
(FPE) for the temporal evolution of their PDF: 



dq d 
dt ~ ~dA 



h 



dg A 
9A ~dA + ^ 



a 2 d ( d g4> dgA a 2 ( d 2 g\q d^ug^q , ^9jfl 



2 dcf> V d(f> dA J 2 I dA 2 dAd<j> dc/) 2 I ' 

(12) 

We are interested in the evolution of A (and h, g A , g<f>, and q are all 2-7r-periodic in <f>), so we integrate both sides of 
the FPE (H2J) from </> = to (f> = 2tt to obtain 

-ot=-BA [h{m + Ta^{J 9Aqd V'YdAj \ 9A ^A +g m) qd ^ (13) 

where Q(A,t) = q(A,<t>,t)d(f>. 

To make further progress, we use an averaging assumption, similar to that in Ref. 

HE!, that thc pDF Q(A,t) 

evolves slowly compared to the time scale of oscillations. Because the mean length of a period is T — 2it/ujo, we 
replace the time-dependent terms in equation (|13j) by their temporal averages over one period of length T. This yields 
the following polynomials: 

1 r T 1 00 

Pi(A) = - g 2 A dt = -(l-A) 2 J2A n \c n \ 2 , 

1 J( > 1 n=l 

P*(A) = ±l gjjfadt+^j^ 9(f! ^±dt = l(l-A)J2A n - 1 (n-A)\c n \\ 



It follows from the forms of the polynomials Pi (A) and P2(A) that that the precise form of the function Z has no effect 
on synchrony for a fixed value of \c n \ 2 , where n is the nth Fourier mode of Z. (This contrasts with the coupling term 
am(9j — 6i), where changing the sinusoidal coupling to other forms of coupling — including cosinusoidal coupling — 
changes the system's qualitative dynamics.) 
The averaged equation for Q is 



2 dA 2 



[P 1 (A)Q(A,t)} . (14) 
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We let te? = to seek a steady-state PDF Qoo(t). Integrating once with respect to A yields 



/*(A) + -P 2 {A) 



Qc 



TdA 



[W)Q°c] , 



(15) 



which can be solved to give the expression in the main text: 



Qoo{A) 



C 



exp 



[2h(A)+a 2 P 2 {A)]dA \ 
° 2 Pi{A) ) 



(16) 



Analysis of Antiferromagnetically Coupled Oscillators in the Absence of Noise 

In the antiferromagnetically coupled system, oscillators are positively coupled to those with the same parity and 
negatively coupled to those with the opposite parity. 

By considering odd oscillators and even oscillators to form separate groups of oscillators, we define complex order 
parameters 



N/2 



N/2 



N ^ 

3 



exp(i02j-i) arid 



N 



(17) 



for the odd and even oscillators, respectively (cf. Ref. [2lj). The antiferromagnetic order parameter defined in the 
main text is r a f = (r c + r e )/2. As with the usual global, equally weighted, sinusoidal coupling in the Kuramoto model 
(which we call ferromagnetic coupling), we let the number of oscillators N — > oo and examine continuum oscillator 
densities f 0}e (uj,6,t). Following the analysis for the Kuramoto model in Ref. [5j, the continuity equations for the 
densities of the oscillators take the form of a pair of nonlinear FPEs: 



dfo,, 
dt 



d_ 

do 



K K 
Y r o,e sin(0) - -^ r e,o sin(0) ) f a>e 



0. 



(18) 



One can then apply Kuramoto's original analysis [3] to this system. Alternatively, one can proceed as in the ferro- 
magnetic case and apply the OA ansatz separately to each family of oscillators. One thereby obtains the coupled 
ordinary differential equations (ODEs) 



dr 0: , 
dt 



-7>"o,i 



r e,o) r o,e( r o,e r e,o)] 



(19) 



Taking the sum and difference of the two equations in (|19l) yields 



d{r e 



dt 

d(r e + r, 



-7 + iuj + 



K 



K 

1 \ , I 2 * , 2 * , 2 * 2 *\ 

[r e ~ r Q ) + — (-r e r e + r e r Q + r r - r r e ) , 



- = -l{r e + r ) + ^r (-r 2 e r* e + r 2 e r* - r 2 r* + r 2 r* e ) 



dt ' 1 4 

In the case with ferromagnetic coupling, we let r = ^/Aexp(uiot 



(20) 



If one were to proceed analogously in 



antiferromagnetic coupling and define r , e = y/ A Q _ e exp(wo£ + 4>o.e), one would obtain four coupled SDEs for A 0iC 
and 4> ,e, and it is then difficult to make analytical progress. However, we seek to quantify the aggregate level of 
synchrony only in the absence of noise. In this case, steady states satisfy A e = A Q and 4> e — —<fi — tp/2, where ip is 
the phase difference between the two groups. (We cannot use this method in the presence of noise, as noise breaks 
the symmetry.) 



Equations (f20|) then simplify to 

' ip\ dtp 



dA . ftp 
—r sin — 
dt \2 

dA ftp 
—— cos — 
dt \ 2 



A cos 
A sin 



2 J dt 

tp\ dip 

2) H 



-27 A sin 
-27 A cos 



KA sin 



-KA 2 



4> 



-KA 2 



. (ty 
sin — 

V 2 



cos 



Sip 
~~2 



cos 



(21) 
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This, in turn, yields 



dA 
~dt 
dip 
~dt 



-27 A + KA(1 - A) sin 2 
-xK{l + A) sinip. 




(22) 



(23) 



By setting ^ — ^ = 0, we seek equilibria of the system. When K sin 2 (ip/2) < 27, there is an unstable equilibrium 
at {A, ip) = (0,0) and a stable equilibrium at (A,ip) — (0,7r). When i^sin 2 ('0/2) > 27, this equilibrium point is 
unstable. Additionally, there is an unstable equilibrium at (A,ij)) = (l — j^ sin a^y 2 ) > ®) an< ^ a s ^ a ble equilibrium at 
(A,ip) = fl — Ksi^(ip/2) ' 7r ) • ^ n P rac tice, this implies that ip(t) — > it, so the threshold for observing synchrony is 

K c = 27 (just as in the Kuramoto model). Similarly, the antiferromagnetic order parameter |r a f| = Vj4sin(^/2) has 
a stable steady state at min {0, y/l - K c /K} ^which has the same dependence on K as the Kuramoto order parameter 
does in the traditional Kuramoto model [J, |5| . We plot the antiferromagnetic order parameter versus the coupling 
strength K in Fig. [5] and obtain excellent agreement with direct numerical simulations of the coupled oscillator system. 
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FIG. 5: Antiferromagnetic order parameter |r a f | versus coupled strength K for width parameter 7 = 1 in the absence of noise 
(i.e., for £7 = 0). The solid curve is the analytical steady state, and circles are from direct numerical simulations of the ODEs 
for an ensemble of N = 500 oscillators. 



